Nuclear stopping and rapidity loss in Au+Au collisions at v / sattv=62.4 GeV 
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Abstract 

Transverse momentum spectra of protons and anti- protons measured in the rapidity range < y < 3.1 from 
0-10% central Au+Au collisions at ^Js^n = 62.4 GeV are presented. The rapidity densities, dN/dy, of protons, 
anti-protons and net-protons (N p — Np) have been deduced from the spectra over a rapidity range wide enough 
to observe the expected maximum net-baryon density. From mid-rapidity to y = 1 the net-proton yield is 
roughly constant (dN/dy ~ 10), but rises to dN/dy ~ 25 at 2.3 < y < 3.1. The mean rapidity loss is 2.01 ±0.16 
units from beam rapidity. The measured rapidity distributions are compared to model predictions. Systematics 
of net-baryon distributions and rapidity loss vs. collision energy are discussed. 
PACS numbers: 25.75 Dw. 



In collisions between gold nuclei at the top en- 
ergy (y/SNN = 200 GeV) of the Relativistic Heavy 
Ion Collider, RHIC, there is strong evidence of a state 
of matter characterized by partonic (quark and gluon) 
degrees of freedom and with properties similar to that 
of a nearly perfect liquid P, 0, H, 0, @. The par- 
tons are produced copiously during the initial stages 
of the collisions and subsequently hadronize into the 
roughly 7000 0] particles produced in central colli- 
sions. The energy required for producing these parti- 
cles comes from the kinetic energy lost by the baryons 
in the colliding nuclei. Since E — mx cosYiy, where 
hit = \f TnP + Pj, is the transverse mass (pr is the 
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transverse momentum) , and y is the rapidity, this en- 
ergy loss is manifest as a loss in mean rapidity of these 
baryons. 

The net-baryon yield can be estimated from the 
net-proton yield, which is taken as the difference be- 
tween the measured yields of protons and anti- protons. 
The rapidity distribution of the net-protons after the 
collision then not only determines the energy available 
for particle production, but also yields information on 
the stopping of the ions due to their mutual interac- 
tion. 

By having a measurement at an energy between 
the lower energy AGS (^/snn = 5 GeV) and SPS 
{^snn — 17 GeV) data and the highest energy RHIC 
results (-^/ijviv = 200 GeV) 



the development 
of rapidity and energy loss can be studied in greater 
detail. In this Letter, we present the first measure- 
ments of the rapidity loss of Au ions after central col- 
lisions at ^/snn — 62.4 GeV. The experimental ar- 
rangement of the BRAHMS detector at RHIC makes 
it possible to measure the distribution of net-protons 
over a rapidity interval y — to y — 3.1. This rapidity 
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range is wide enough to include the maximum proton 
rapidity density, in contrast to the situation at the 
RHIC top energy, where the beam rapidity is higher 
(yb = 5.4 compared to 4.2 at 62 GeV) and where the 
acceptance of existing experiments does not include 
this peak. The situation at the lower energy there- 
fore makes it possible to better determine the rapidity 
density distribution. Together with similar informa- 
tion from experiments at lower energies, and at the 
RHIC top energy, we conclude that the mean rapidity 
loss of ultra-relativistic heavy ions exhibits a slowly 
varying behavior as a function of beam energy from 
SPS energies on upwards. 

The BRAHMS detector consists of two magnetic 
spectrometers: the Mid Rapidity Spectrometer (MRS) 
able to cover polar angles (measured with respect to 
the beam direction) 30° < 8 < 90° and the Forward 
Spectrometer (FS) able to cover 2.3° < 6 < 15°. Each 
spectrometer determines trajectories and momenta of 
charged hadrons. Two Time Projection Chambers 
(TPCs) are utilized in the MRS and two TPCs and 
three Drift Chambers (DCs) in the FS. Together, they 
measure protons and anti-protons in the range — 0.1 < 
y < 3.5. Collision centrality is determined using a sil- 
icon and plastic tile multiplicity array located around 
the nominal intersection point (NIP) [6]. For this 
analysis, the centrality class — 10% was selected, 
corresponding to a calculated number of participant 
nucleons of N part = 314±8. The vertex position is de- 
termined with an accuracy of ~ 1 cm using two arrays 
of Cherenkov counters positioned on either side of the 
NIP [lflj. We have selected for analysis tracks in the 
MRS (FS) with vertices within ±15(20)cm from the 
NIP. 
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Figure 1: Spectra of identified protons and anti-protons for 
y ~ 0, y ~ 0.65, y ~ 2.3 and y ~ 3 respectively. The solid 
drawn lines are the fit functions used to determine the yield. 
Vertical bars show statistical errors only. 

Particle identification (PID) is done in the MRS 
via time of flight (TOF) measurements, enabling clean 
identification of protons and anti-protons in the mo- 
mentum range 0.4 GeV/c < p < 3 GeV/c. In the 



FS, two PID detectors are used: A TOF system and 
a Ring Imaging Cherenkov (RICH) detector. In the 
RICH, only (anti-)protons with p > 15 GeV/c will 
create resolvable rings. Lower momentum (anti)protons 
have no associated Cherenkov radiation. In the mo- 
mentum range 3 GeV/c < p < 7.5 GeV/c we apply a 
2 sigma cut about the (anti-) proton peak in the cal- 
culated mass spectrum based on the TOF measure- 
ment and require that the particle not be identified in 
the RICH as a pion. In the range 12 GeV/c < p < 
20 GeV/c the RICH is used for PID. For 12 GeV/c < 
p < 16 GeV/c the (anti)protons are those particles 
with either no associated ring or a ring with a small 
radius. Since the RICH is 97% efficient for particles 
above threshold, up to 3% of pions and kaons in this 
momentum range will be mistakenly identified as pro- 
tons, an effect for which we correct the data. 
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Figure 2: Top: rapidity densities of protons and anti-protons. 
Bottom: rapidity densities of net protons. The statistical er- 
rors are smaller than the marker sizes. The verticalal brackets 
shows the estimated systematic errors due the extrapolation of 
yields, and the horizontal bars indicate the width of the ra- 
pidity interval for each data point. The histograms are from 
HIJING/BB(v2.1). 

The differential invariant yields, 2 ^ pT J^J^ T , of pro- 
tons and anti-protons have been corrected for geomet- 
rical acceptance and detector efficiency. The accep- 
tance correction is, due to the small solid angle of the 
BRAHMS spectrometers, the largest correction and is 
obtained from a purely geometrical simulation of each 
combination of angle and magnetic field used in the 
experiment. The PID efficiency of the TOF walls is 
found to be 93-98%. The total tracking efficiency in 
the MRS and FS is » 90% and » 80%, respectively. 
The data have also been corrected for absorption, en- 
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ergy loss and multiple scattering estimated through 
simulations using GEANT The corrections for 

these physics effects amount to less than 20% at the 
lowest pt and less than 15% at the highest pr for the 
FS, whereas in the MRS the correction at the highest 
p T is 8%. 

Fig. [T] shows the invariant spectra of protons and 
anti-protons for y ~ 0, y ~ 0.65, y ~ 2.3 and y ~ 3 
as a function of pt- The integrated yields have been 
obtained by fitting these spectra with Gaussian func- 
tions (/(pr) of e:rp(— p T /2cr 2 )) and integrating the 
fit functions in the range < pr < oo. From the 
yields we obtain the rapidity densities ^p. The inte- 
grated yield under the data points compared to the 
extrapolated yield is - 80 - 90% in the MRS and ~ 
40 — 50% in the FS. The mean pt and tot determined 
from the distributions varies from (jpr) ~ 1 GeV/c 
((to t ) - 1.4 GeV) in the MRS to (p T ) ~ 0.7 GeV/c 
((to t ) - 1.2 GeV)in the FS. 

The top panel of Fig.[2]shows the dN/dy of protons 
and anti-protons. The lower panel shows the dN/dy 
of the net-protons, defined here as the difference of 
the protons and anti-proton densities. The system- 
atic errors due to the limited pr coverage have been 
determined by varying the fit range, the bin width, 
the choice of fit function and the rapidity cuts. These 
uncertainties together with those for the corrections 
results in an estimate of the total systematic error of 
roughly 10%. 

The predictions of the HIJING/BB(v2.1) event gen- 
erator [l2j for the proton, anti-proton and net-proton 
rapidity distributions are shown in Fig. [2] by the his- 
tograms. The HIJING distributions include protons 
and anti-protons from hyperon decays that would be 
measured as primary protons (e.g. 0.53 ± 0.05 (anti- 
protons for each (anti-)lambda, see The HIJING 
calculation predicts a smaller net-proton yield at mid- 
rapidity, and shows a peak slightly forward of that in- 
dicated by our results; thus, the mean rapidity loss is 
smaller than estimated from our data. 

The estimate for net-baryons depends on the rel- 
ative acceptance for direct and decay protons in the 
spectrometers and the ratio of hyperons to protons. 
The acceptances are calculated using GEANT while 
the hyperon:proton ratios have been calculated us- 
ing the THERMUS [H model. The THERMUS re- 
sults are based on fits to BRAHMS 7r ± ,i ; r ± ,p and p 
data 0. The n/p ratio is obtained from HIJING. 
The necessary conversion from the net-protons, Np, 
to the net-baryons, Nb, is then done (as in Ref. [9j) 
using Nb = (2 ± 0.1) ■ N p at mid-rapidity and Nb = 
(2.1 ± 0.1) • N p at forward rapidities (the larger cor- 
rection at forward rapidity is due to a small increase 
in the n/p ratio at forward rapidities). 

We quantify the mean rapidity [l5| and energy 



loss using the following integrals: 
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Here yb and E are the rapidity and energy of the 
incoming beams and dN B -B/dy is the net-baryon ra- 
pidity density. Since we only have data at four rapidi- 
ties, y = 0, 0.65, 2.3 and 3, we have to both interpolate 
and extrapolate our data to evaluate the integrals. We 
found that (tot) drops linearly with rapidity. For the 
baryon yield we fitted our data to a 3 rd order polyno- 
mial in y 2 (analogous to [9() subject to the constraint 
that Np ar t = 314 and that the yield at yb is 0. The 
result of this fit together with our net-baryon distri- 
bution arc shown in the inset of Fig. [3J 

Using this method we obtain the average rapidity 
and energy loss per participant baryon for Au+Au 
collisions at ^snn = 62.4 GeV: 

Sy = 2.01 ± 0.16, SE = 22 ±1 GeV. 

Since the energy loss per pair of participant nucleons 
is 2SE, approximately 70% of the initial beam energy 
is available for particle production after the collision. 
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Figure 3: Rapidity losses from AGS, SPS and RHIC as a func- 
tion of beam rapidity. The solid line is a fit to SPS and RHIC 
data, and the band is the statistical uncertainty of this fit. The 
dashed line is a linear fit to AGS and SPS data from [15)1 . 

The rapidity losses at AGS 0,0 E3|, SPS d and 
RHIC y/Fi^ = 200 GeV [§] together with the present 
result are summarized in Fig. [3] The figure shows that 
rapidity loss increases rapidly with beam energy from 
AGS to SPS but much more slowly from SPS energy to 
RHIC energies. The dotted line in Fig.[3]is taken from 
(la l where there was found to be a linear scaling over a 
wide range of energies up to the SPS top energy. This 
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scaling was found to be broken at -^snn — 200 GeV in 
Q and we conclude that the proposed linear scaling 
breaks between ^snn — 17 GeV and ^/snn — 62.4 
GeV. The solid drawn line is a linear fit to SPS and 
RHIC data and allows one to extrapolate to the LHC 
regime (j/f, = 8.7). The grey band gives the statistical 
uncertainty of this extrapolation, but is only useful 
to the extent that the underlying physics is the same 
from RHIC to LHC energies. 
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dN/dy' for net-baryons normalized to the number of 
participant pairs. For NA49 data ;8f] the number of 
net-baryons is 352 ± 12 for the interval \y cm \ < 2.5. 
Extrapolating to j/f, = 2.9 we estimate the number of 
participants to be 390 ± 20. 

Fig. [4] (top left) shows net-brayon rapidity densi- 
ties \ / (N part /2)dN B _jj / dy 1 for different beam ener- 
gies. The y/SNN = 17 GeV and ^/snn — 62.4 GeV 
data differ significantly at y' ~ — 3 but approach each 
other towards yb, coinciding at y 1 ~ — 1, reflecting 
the much larger contribution of 'target' baryons at 
\J s nn = 17 GeV. This is reminiscent of limiting frag- 
mentation and indeed is the reason for choosing the 
variable y' . 
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Figure 4: The top left panel shows (1/N part /2)dN B _g /dy' 
where y' = y — yb for SPS and RHIC energies. No scaling 
is observed. The three other panels show data from SPS and 
RHIC drawn together with the 'target' net-baryon contribution 
from 11911 ■ The triangles represent data from NA49, the squares 
BRAHMS data from 62 GeV and the circles BRAHMS 200 GeV 
data. 

Applying the conversion from net-protons to net- 
baryons and fitting the result, we estimate that for col- 
lisions with N part ~ 314, there are ~ 240 net-baryons 
in the acceptance covered by BRAHMS, —3.1 < y < 
3.1 (because of symmetry). To estimate the minimum 
and maximum possible stopping at v/sjvjv = 62.4 
GeV the remaining baryons can be placed at y — 2.7 
or y = y b . These values are indicated by the hor- 
izontal bars in Fig. [3J The horizontal bars for the 
y/SNN — 200 GeV are obtained in the same way, and 
similarly give the absolute limits on the amount of 
stopping. If instead of relying on these limits, we use 
the uncertainty on the fits to the yields to perform 
the extrapolation to LHC, we then find that the ra- 
pidity loss at LHC would be expected to lie between 
2.1 < Sy < 2.4. 

The slow increase of the rapidity loss from the top 
SPS energy to the top RHIC energy indicates that 
baryon transport does not depend strongly on energy 
at high energies when observed in the rapidity frame 
of the beam, y' — ycM — Vb- We therefore compare 



Figure 5: Projectile net-baryon rapidity density 

(l/N part /2)dN^ r ° J £ ctile /dy' from SPS and RHIC after 
subtraction of the target net-baryon contribution (see Fig.|lJ. 

This direct comparison of data is complicated by 
the aforementioned target contribution. For ^snn = 
17 GeV at y' = 3 the target contribution is half the 
net-baryons while for ^snn — 62.4 GeV and ^snn = 
200 GeV the contribution is significantly less. To com- 
pare the net-baryons from the projectile only, the tar- 
get contribution must be subtracted. At ycM = 
the target nucleons comprise half of the net-baryons 
due to beam-target symmetry, and at y^ the target 
contribution is expected to be negligible. Between 
these extremes we consider two different rapidity de- 
pendences for the correction: (1) a simple exponential 
form exp(— y') [18 1 and (2) a gluon junction motivated 
form exp(— y'/2) 19(. These two functions give limits 
for the target net-baryon tail and provide the bounds 
for the grey bands in Fig. |U We subtract the av- 
erage of these two functions, with half the difference 
between them taken to be the systematic uncertainty 
on the function, from the data. A consequence of the 
forward-backward (beam-projectile) symmetry is that 
there is no error on the beam contribution at mid- 
rapidity. As one approaches the projectile rapidity, 
the uncertainty grows, but both the absolute and rel- 
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ative contribution of the target net-baryons decreases 
and the correction becomes smaller. After subtract- 
ing the target net-baryon contribution shown in Fig. [H 
we obtain the projectile net-baryon distributions pre- 
sented in Fig. [5l One notes a remarkable similarity 
among the data sets for large rapidity losses y' < —3, 
while the 200 GeV data may begin to diverge from 
the lower energy data for y' > —3. This suggests that 
the rapidity loss obtained from Eq. 1, which does not 
distinguish the target and projectile contributions to 
the net-baryon yield, while experimentally practical, 
is slightly misleading, and that the increase in projec- 
tile rapidity loss from SPS to RHIC maximum energy 
might be smaller than what is implied in Fig. [3] 

Finally, we note that the present results for the ra- 
pidity loss of baryons (protons) for Au+Au collisions 
at ^/snn — 62.4 GeV, together with similar data at 
lower and higher energy permit us to more accurately 
determine the onset of a collision regime where the ab- 
solute rapidity loss appears to vary slowly with beam 
rapidity with the data exhibiting a near saturation of 
the rapidity loss with increasing collision energy. Re- 
cently, it has been argued [20( that such a scaling may 
result from gluon saturation of the collision. These 
results suggest that the rapidity loss at LHC will be 
approximately Sy = 2.3. Although the baryon rapid- 
ity loss relative to the beam rapidity decreases, the 
overall energy available for particle production still 
increases with increasing beam energy. 
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